1/1 


MICROCOPY  RESOLUTION  TEST  CHART 

■  JRMO  u- 


I 


NAVAL  POSTGRADUATE  SCHOOL 
Monterey ,  California 


c 


DTIC 

^ELECTE 
MAY  3  1 1988 


I 


CALIBRATION  OF  A  TRIPLE  WIRE  PROBE 
FOR 

TURBULENCE  MEASURMENTS 


Gary  J.  Selman 


December  1987 


Thesis  Advisor 


J.  Val  Healey 


Approved  for  public  release;  distribution  is  unlimited. 


8  8  5  27  081' 


J* -J..'  Jii'iH 


URITY  CLASSIFICATION  OF  THIS  PAGE 


la  REPORT  SECURITY  CLASSIFICATION 

Unclassified 


2a.  security  classification  authority 


2b  OECLASSIFICATION/ DOWNGRADING  SCHEDULE 


REPORT  DOCUMENTATION  PAGE 


lb  RESTRICTIVE  MARKINGS 


3  DISTRIBUTION /AVAILABILITY  OF  REPORT 

Approved  for  public  release,  distribution 
unlimited 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

S  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

6a  NAME  OF  PERFORMING  ORGANIZATION 

6b  OFFICE  SYMBOL 
(If  applicable) 

7a.  NAME  OF  MONITORING  ORGANIZATION 

Naval  Postgraduate  School 

CODE  67 

Naval  Postgraduate  School 

6c.  ADDRESS  (City.  Stale,  and  ZIP  Code) 

7b.  ADDRESS  (Crty,  State,  and  ZIP  Code) 

Monterey  CA,  93943  -  5100 

Monterey  CA,  93940-5100 

8a  NAME  OF  FUNDING  /SPONSORING 
ORGANIZATION 

8b  OFFICE  SYMBOL 
(If  applicable) 

9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

8c.  ADDRESS  (City.  State,  and  ZIP  Code) 

<0.  SOURCE  OF  FUNDING  NUMBERS 

PROGRAM 
ELEMENT  NO. 


PROJECT 

NO 


WORK  UNIT 
ACCESSION  NO. 


II  TITLE  (Include  Security  Classification) 

Calibration  of  a  Triple  wire  Probe  for  Turbulence  Measurements 


12  PERSONAL  AUTMOR(S) 


13a  TYPE  OF  REPORT 

Masters  Thesis 


Gary  J.  Selman 


13b  TIME  COVERED 
FROM _ TO 


14  DATE  OF  REPORT  (Year.  Month.  Day)  IIS  PAGE  COUNT 
1987,  December  1987  I  44 


17 

COSATI  CODES 

FIELD 

GROUP 

SUB-GROUP 

18.  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 


Hot-Wire i  Triple  Wire  Probe 


19.  ABSTRACT  ( Continue  on  reverse  if  necessary  and  identify  by  block  number) 


A  method  of  calibration  for  a  triple  wire  probe  has  been  implimented 
for  use  in  the  flow  visualization  tunnel  at  the  •  Naval  Postgraduate 
School.  The  system  was  designed  to  allow  flow  measurements  in  highly 
turbulent  flows  and  under  conditions  of  relatively  unknown  flow 
directions.  Preliminary  evaluation  of  the  method  indicates  that 
accuracies  in  the  range  of  1-4%  can  be  expected.  Although  attainment  of 
these  levels  of  accuracy  have  not  yet  been  achieved,  it  is  felt  that 
they  can  be  realized. 


20  DISTRIBUTION  /AVAILABILITY  OF  ABSTRACT  21.  ABSTRACT  SECURITY  CLASSIFICATION 

□  UNCLASSIFIEO/UNLIMITEO  □  SAME  AS  RPT  □  DTIC  USERS 


22».  NAME  OF  RESPONSIBLE  INDIVIDUAL  1 1 2b.  TELEPHONE  (Include  Area  Code) 

.7.  Val  Mealy _  I  (408)  646-2491 


DO  FORM  1 473, 84  mar  83  APR  edition  may  be  used  until  exhausted.  un  ...tv  r,  .c( 

All  Other  editions  «re  obsolete  - SVA>? 


\  PTHHHif  omesi  I  ft  HI  M. 


1 


Approved  for  public  release;  distribution  is  unlimited. 


Calibration  of  a  Triple  Wire  Probe 
for 

Turbulence  Measurments 


by 


Gary  J.  Selman 
Lcdr,  United  States  Navy 
B.S.,  The  Citadel,  1977 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  AERONAUTICAL  ENGINEERING 


from  the 


NAVAL  POSTGRADUATE  SCHOOL 
December  1987 


ABSTRACT 


A  method  of  calibration  for  a  triple  wire  probe  has  been 
implemented  for  use  in  the  flow  visualization  tunnel  at  the 
Naval  Postgraduate  School.  The  system  was  designed  to  allow 
flow  measurements  in  highly  turbulent  flows  and  under 
conditions  of  relatively  unknown  flow  directions. 
Preliminary  evaluation  of  the  method  indicates  that 
accuracies  in  the  range  of  1  -  4%  can  be  expected.  Although 
attainment  of  these  levels  of  accuracy  have  not  yet  been 
achieved,  it  is  felt  that  they  can  be  realized. 


j^Acoesslon  For 

~Nt7s  ~GRA*I 
DTIC  TAB 
Unannounced 
Just Ifl oat  Ion. 


3 


i  By - 

Plat r i but  1 on/ _ 

Availability  Codas 


(Diet 


Avail  and/or 
Special 


□  □ 


TABLE  OF  CONTENTS 


I.  INTRODUCTION  .  7 

I.  EQUIPMENT .  9 

A.  PROBE .  9 

B.  DATA  COLLECTION  EQUIPMENT  .  9 

C.  CALIBRATION  EQUIPMENT  .  10 

II.  PROBE  CALIBRATION  .  15 

A.  GEOMETRIC  DEFINITIONS  .  16 

B.  CALIBRATION  THEORY  .  16 

C.  DATA  REDUCTION . 20 

D.  MODIFICATIONS  TO  BASIC  PROCEDURE  ....  22 

III.  SOFTWARE  GUIDE . 23 

A.  CALIBRATION  PROCEDURES  .  23 

B.  DATA  COLLECTION . 25 

C.  DATA  REDUCTION . 26 

D.  GENERAL  NOTES  ON  THE  SOFTWARE . 26 

IV.  CONCLUSIONS . 28 

V.  RECOMMENDATIONS . 29 

APPENDIX . 31 

REFERENCES . 42 

INITIAL  DISTRIBUTION  LIST . 43 


4 


LIST  OF  FIGURES 


1.  Velocity  Calibration  .  12 

* 

2.  Probe  Positioner  .  13 

3.  Velocity  Profile  .  14 

4.  Coordinate  System  .  17 

5.  Velocity  Magnitude  Calibration  .  21 

6.  Wire  Angle  Calibration . 21 


\V  v 


V.v 


ACKNOWLEDGEMENTS 

First  and  foremost,  I  would  like  to  thank  my  loving 
wife,  Lynda,  without  whose  support  the  completion  of  this 
research  would  have  been  impossible.  Her  understanding  of 
all  those  late  nights  and  her  comfort  when  things  went  wrong 
will  always  be  remembered. 

I  would  also  like  to  thank  my  Thesis  Advisor,  Professor 
Val  Healey,  who  provided  me  with  the  opportunity  to  do  this 
research . 

Last  but  not  least,  I  thank  Lcdr.  Tom  Cahill,  who  by  his 
help  in  areas  outside  of  this  research,  provided  me  the  time 
to  complete  this  project. 


I.  INTRODUCTION 


A  method  of  calibration  of  a  three  component  constant 
temperature  anemometer  (CTA)  and  the  reduction  of  data 
obtained  from  the  system  is  discussed.  This  report  is 
intended  to  document  the  system  for  future  reference  in  the 
study  of  flow  around  ship  models  currently  being  conducted 
at  the  Naval  Postgraduate  School  (NPS). 

Recent  work  in  the  Flow  Visualization  Tunnel  at  NPS  has 
successfully  modeled  the  atmospheric  boundary  layer  that 
exists  in  open  ocean  conditions  [Ref.  1].  Extensive  on-  and 
off-body  flow  visualization  studies  have  also  been  conducted 
to  determine  the  flow  field  around  a  generic  destroyer 
model.  Work  currently  being  completed  will  provide  the 
characteristics  of  the  flow  around  an  oscillating  model. 
These  studies  are  being  conducted  as  part  of  a  study  to 
determine  the  feasibility  of  simulating  the  ship/helicopter 
interface  with  sufficient  realism  to  permit  certification  of 
helicopter  operations  on  those  ship  decks  not  currently 


The  next  step  in  this  study  is  to  map  the  flow  field 


around  the  model.  The  three  components  of  mean  velocity, 
flow  directions,  turbulence  intensities,  and  spectrum 
functions  are  desired.  The  above-mentioned  work  has 
demonstrated  the  extremely  turbulent  nature  of  the  flow  as 
well  as  the  locations  of  reverse  flow. 

Historically,  analog  circuits  have  been  used  to  reduce 
wire  voltages  to  velocity  components.  This  has  been  due, in 
part,  to  the  amount  of  computation  time  required  per  data 
point.  Butler  and  Wagner  [Ref.  2]  and  Westphal  and  Metha 
[Ref .  3  ]  have  derived  methods  of  calibration  and 
measurement  that  have  achieved  a  considerable  reduction  in 
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required  computational  and  calibration  time. 


A .  PROBE 


A  commercial  DISA  triaxial  probe,  type  55P91  was  used 
for  measurements.  The  wires  of  the  probe  are  arranged 
parallel  to  an  orthogonal  coordinate  axis  system  that  is 
rotated  with  respect  to  the  probe  axis  system.  The  ends  of 
the  wires  lie  on  the  surface  of  a  sphere  that  has  a  diameter 
of  3  mm,  so  that  the  sensitive  volume  of  the  probe  is 
smaller  than  that  of  the  9  mm  diameter  of  the  head  [Ref.  4, 
p.  17].  The  reference  position  of  the  probe  is  with  itB 
number  three  wire  in  the  X-Y  plane  with  its  shortest  prong 
at  the  smaller  value  of  Y. 

B.  DATA  COLLECTION  EQUIPMENT 

A  three  channel  Dantec  constant  temperature  anemometer 
system  was  used  for  data  collection.  Each  channel  consisted 
of  a  56C01  CTA  with  a  56C17  CTA  Bridge  Module.  The  output  of 
the  CTA  was  further  processed  by  a  56N20  Signal  conditioner 
that  allowed  filtering  of  unwanted  frequency  components. 
Also  available  in  the  system  was  a  56N25  Digital  RMS 


voltmeter  and  a  56N22  Digital  Mean  Value  Unit.  The  latter 
two  units  were  used  to  monitor  performance  of  system 
software  and  hardware  during  flow  measurement. 

The  outputs  from  the  three  CTA  channels  were  connected 
to  a  12  bit  Metrabyte  A/D  converter  where  they  were 
collected  and  stored  for  reduction.  The  A/D  board  uses  a  12 
bit  successive  approximation  converter  with  a  12  microsecond 
conversion  time  yielding  a  maximum  throughput  rate  of  60  KHz 
in  D.M.A  mode  [Ref.  5,  p.  1].  Analog  inputs  are  sampled  in 
sequence  by  means  of  a  multiplexer  which  passes  the  data  to 
a  sample  and  hold  circuit  where  they  are  held  until 
completion  of  the  previous  conversion.  This  causes  a  delay 
of  approximately  one  conversion  time  between  successfully 
scanned  channels.  Due  the  relatively  low  frequency  nature  of 
the  signals  to  be  studied  (below  lKHz),  This  delay  was  not 
considered  significant. 

C.  CALIBRATION  EQUIPMENT 

Calibration  was  conducted  in  a  small  wind  tunnel  with  a 
5"  X  5"  interior  cross  section.  Static  pressure  ports  were 
located  along  one  side  of  the  tunnel  to  provide  test  section 
velocities.  The  static  pressure  was  measured  by  means  of  a 
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micromanometer  with  one  port  connected  to  a  static  pressure 
port  opposite  the  probe  head.  The  other  port  was  open  to 
atmospheric  pressure.  Actual  velocity  was  derived  from  a 
calibration  curve  relating  the  static  pressure  differential 
to  velocities  obtained  from  a  pitot-static  tube  (Fig.  1). 
Due  to  the  small  size  of  the  tunnel,  it  was  undesirable  to 
mount  a  pitot-static  probe  in  the  tunnel  during  calibration. 

The  probe  positioner  was  designed  to  provide  ±  45*  of 
pitch  and  ±  90*  of  yaw.  and  was  mounted  on  top  of  the 
calibration  tunnel.  Figure  2  shows  a  schematic  of  the 
positioner  and  probe  holder.  The  pivot  point  for  pitch  and 
yaw  were  coincident  and  are  located  on  the  center  line  of 
the  tunnel.  Yaw  was  achieved  by  means  of  a  stepper  motor 
under  the  control  of  the  calibration  program.  This  stepper 
motor  allowed  for  yaw  increments  of  15*  per  step.  During 
calibration,  pitch  was  adjusted  manually.  Markings  of  ± 
5,10,20,30,40,  and  45’  were  etched  on  the  positioner  arch. 

The  stem  of  the  probe  holder  was  offset  so  that  the  head 
of  the  probe  was  approximately  1.5"  ahead  of  the  pivot  point 
of  the  system.  This  arrangement  allowed  for  maximum  pitch 
and  yaw  ranges  while  minimising  probe  head  travel.  A  survey 
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III.  PROBE  CALIBRATION 


As  previously  mentioned,  a  commercially  available  probe 
having  three  independent,  mutually  orthogonal  wires  was 
utilized.  A  method  of  calibration  was  sought  that  provided 
rapid  reduction  of  data  and  reduced  calibration  time. 
Several  methods  of  calibration  were  investigated  that 
provided  highly  reliable  measurements  of  mean  and 
fluctuating  velocity  components.  However,  in  flows  of 
relatively  large  turbulence  intensity  and  unknown  flow 
directions , these  methods  tend  produce  large  errors  [Ref.  2, 
p.  1].  As  pointed  out  in  the  introduction,  the  calibration 
process  had  to  provide  accurate  results  when  subjected  to 
this  type  of  flow  condition.  A  search  of  the  current 
literature  indicated  that  the  method  of  Butler  and  Wagner 
[Ref.  2]  was  well  suited  to  the  problem  at  hand.  Not  only 
did  it  extend  the  yaw  and  pitch  envelope  but  it  also  provide 
for  a  relatively  rapid  calibration  process.  This  method  is 
discussed  in  detail  in  Reference  1  and  an  overview  with 
slight  modifications  is  given  here. 
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A.  GEOMETRIC  DEFINITIONS 

Wire,  probe,  and  Laboratory  coordinate  systems  are 
defined  in  Figure  4.  Pitch  (a)  is  defined  as  the  angle 
between  projection  of  the  velocity  vector  onto  the  X-Z  plane 
and  the  X-axis.  Yaw  (p)  is  the  angle  between  the  projection 
of  the  vector  in  the  X-Y  plane  and  the  X  axis.  The  positive 
sense  of  these  two  angles  is  as  shown  in  the  diagram.  The 
angle  0 i  is  the  angle  between  the  velocity  vector  and  wire  i 
where  i  =  1,2,3. 

B.  CALIBRATION  THEORY 

The  underlying  theory  of  the  calibration  procedure  is 
that  if  the  angle  between  any  two  wires  of  the  probe  and  the 
velocity  vector  is  known,  the  pitch  and  yaw  of  the  velocity 
vector  in  laboratory  coordinates  can  be  derived.  Using  this 
method,  there  are  two  possible  directions  of  the  velocity 
vector.  The  angle  between  the  third  wire  and  the  velocity 
vector  is  used  to  resolve  the  proper  direction. 

If  the  components  of  a  vector  are  known  in  the  wire 
coordinates,  then  the  components  of  that  vector  in  the 
laboratory  coordinate  system  can  be  found  by  the  following 
equation : 
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Figure  4.  Coordinate  System 
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Xi,Yi,Zi  *  Components  in  laboratory  system 

a,p  *  pitch  and  yaw  of  probe  axis  Eqn.  1 

XwtYwtZw  s  Components  in  wire  system 

Likewise,  the  velocity  vector  can  also  be  described  in 

terms  of  the  laboratory  coordinate  system: 

Q  =  -Q( cosa ) ( cosp ) i  -  Q(sinp)j  -  Q( sina ) ( cosp )k  Eqn.  2 

i,j,k  *  unit  vectors  ;  Q  »  Velocity  Magnitude 

An  expression  for  the  angle  between  the  velocity  vector 

and  any  wire  unit  vector  can  be  found  by  taking  the  dot 

product  of  the  two: 

cos  0i  =  -/3/3(cosa) (cosp)  -  /6/6(sinp)  -  /2/2 ( sina )( cosp ) 

cos  02  =  -V'3/3  ( cosa )( cosp  )  -  /6/6(sinp)  +  /2/2  (  sina )  ( cosp ) 

cos  ©j  =  -/3/3 (cosa ) (cosp )  +  /6/3(sina)  Eqn.  3 

The  following  relationship  was  used  to  relate  the  wire 

voltages  to  the  effective  cooling  velocity: 

Q*  f  t  t  -  K i  +  KiEi  +  KsE*  Eqn.  4 

The  effective  cooling  velocity  was  taken  to  be  equal  to 

the  total  velocity  when  the  wire  was  oriented  normal  to  the 


calibration  flow.  The  coefficients  of  Equation  3  for  each 
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wire  were  determined  by  varying  the  tunnel  velocity  with 


each  sensor  in  this  orientation  and  then  performing  a  least 


squares  fit  on  the  data, 


Probe  angular  response  was  calibrated  by  varying  the 


pitch  and  yaw  of  the  probe  with  respect  to  the  centerline  of 


the  tunnel.  Yaw  and  pitch  angles  as  well  as  the  three  wire 


several  different  tunnel  velocities.  The  range  of  yaw  and 


pitch  were  as  follows 


-30  <  p  <  30  deg 
-40  <  a  <  40  deg 


The  following  correlating  parameters  were  then  used  to 


obtain  accurate  data  fits  for  the  velocity  magnitude  and 


sensor  angles  [Ref.  2,  p.  728J  : 


X  *  V  (Ei3  +  Ei*  +  Ei*  ) 


Velocity  Magnitude  Parameter 


Y  *  E  i  /  /  (E,*  +  Ej*  ) 


Pitch  Angle  Parameter 


Z  s  Ei  /  E: 


Yaw  Angle  Parameter 


The  Calibration  data  for  velocity  magnitude  were  then 


least  squares  fitted  to  the  following  equation: 


Q  =  Ao  +  AxX  ♦  AiX*  +  AiY  ♦  A«Y* 

+  A sZ  +  A*Z*  ♦  A tXY  ♦  A.XZ  +  A*YZ 


Eqn.  5 


Calibration  for  sensor  angle  were  then  fit  according  to 


the  following  : 
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0i  =  Bo  +  BiQ  +  BaQ2  +  BaFi  ♦  BoFi* 

+  BsY  +  B*Y2  +  B7Z  +  BsZ2+  B.QFt 
+  B  i  oQY  +  BnQZ  +  BiiFiY  +  BiaFiZ  ♦  Bi4YZ 


Eqn .  6 


The  least  squares  curve  fitting  of  the  above  two 


equations  was  accomplished  utilizing  MATLAB,  a  commercially 


available  program  by  Mathworks.  The  use  of  this  program 


greatly  simplified  the  programming  effort  as  well  as 


providing  extremely  rapid  data  reduction.  Calibration  data 


for  one  wire  is  shown  in  Figures  5  and  6. 


C.  DATA  REDUCTION 


Once  the  simultaneous  sensor  output  voltages  have  been 


collected  from  the  wires,  the  three  parameters  given  in 


Section  B  are  computed  and  substituted  into  Equation  5.  The 


velocity  magnitude  is  then  utilized  in  Equation  6  to  obtain 


the  wire  angle  for  each  of  the  three  wires. 


As  previously  noted,  there  are  two  possible  solutions 


for  the  pitch  and  yaw  angle  of  the  velocity  vector  given  any 


two  of  the  wire  angles.  The  details  of  the  procedure  to 


resolve  the  proper  solution  is  given  in  the  Appendix  of 


Reference  2. 
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D.  MODIFICATIONS  TO  BASIC  PROCEDURE 

It  should  be  noted  that  the  definitions  of  a  and  b  have 
been  reversed  from  those  given  in  Reference  2.  The  origin  of 
the  pitch  angle  has  also  been  modified.  This  was  done  to 
bring  these  definitions  more  into  line  with  current 
aeronautical  usage.  The  appropriate  changes  to  Reference  2 
equations  were  made  to  account  for  this  change. 

Reference  6  discusses  two  methods  by  which  the  accuracy 
of  hot  wire  data  obtained  using  this  particular  probe  may  be 
improved : 

(1)  Adjustment  of  the  overheat  control  of  two  of  the 
wires  to  match  the  calibration  of  the  third. 

(2)  Compensation  for  temperature  drift  away  from 
calibration  conditions  during  measurement  by 
adjustment  of  overheat  controls  to  produce  the 
same  output  voltages  at  a  given  tunnel  velocity 
as  was  obtained  during  calibration. 

Method  1  was  used  during  testing.  Due  to  the  temperature 
stability  of  the  laboratory  facility,  the  second  method  has 
not  been  utilized. 
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IV.  SOFTWARB  GUIDE 


A.  CALIBRATION  PROCEDURES 

Probe  calibration  consists  of  two  separate  procedures, 
the  calibration  of  each  individual  wire  and  velocity 
magnitude/yaw-pitch  calibration.  The  following  steps  will 
ensure  that  the  necessary  procedures  will  be  carried  out  in 
the  proper  order. 

Step  1 . 

Ensure  that  no  memory  resident  programs  are  currently 
running  as  this  will  interfere  with  the  proper  operation  of 
some  of  the  data  fitting  routines.  Start  the  calibration 
routine  by  typing  "CALIBRAT"  at  the  DOS  prompt.  This  will 
provide  the  user  with  a  menu  of  choices. 

Step  2. 

Each  individual  wire  must  first  be  calibrated  to  obtain 
the  coefficients  of  Equation  4.  Selection  of  option  1,  " 
King's  Law  Calibrate",  will  start  this  portion  of  the 
calibration  process.  Each  wire  in  turn  must  be  placed 
perpendicular  to  the  calibration  flow  during  this  procedure. 
This  involves  a  rotation  of  the  probe  about  its  X-axis  and  a 


yaw  about  the  Z-axis.  The  X-axis  rotation  is  done  manually 
and  is  accomplished  by  loosening  the  set  screw  that  holds 
the  probe  in  the  probe  holder.  Proper  positioning  is 
accomplished  by  locating  one  of  the  three  pattern  of  dots 
that  indicate  the  individual  wires  on  top  of  the  probe  axis. 
The  software  calibration  routine  provides  on  screen 
documentation  of  the  procedure.  After  positioning  the  probe 
according  to  the  directions,  the  software  yaws  the  probe  so 
as  to  bring  the  wire  under  calibration  perpendicular  to  the 
flow.  After  the  calibration  data  has  been  collected,  the 
procedure  is  repeated  until  all  wires  have  been  calibrated. 

The  software  will  generate  a  file  called  "VELCAL.M" 
which  contains  the  data  for  each  wire.  Data  is  stored  with 
the  voltage  in  Column  1  and  the  calibration  velocity  in 
Column  2.  This  data  is  then  fit  by  a  MATLAB  file  called 
"DATAFIT".  The  results  are  stored  in  a  file  titled  "KPARAM" 
which  is  later  used  by  the  data  reduction  program  as  well  as 
the  MATLAB  file  that  performs  the  curve  fitting  on  the 
velocity  and  wire  angle  data.  Entry  into  the  matlab  routines 
is  automatic  and  no  user  interaction  is  required. 


Data  for  velocity  magnitude  and  wire  angle  calibration 
is  collected  simultaneously.  A  minimum  of  15  data  points  at 
various  yaw, pitch  and  velocity  settings  are  required  to 
complete  the  process.  After  selection  of  option  2, 
"Yaw/Pitch  Calibrate",  answer  the  questions  as  they  are 
asked.  Probe  pitch  angle  must  be  set  manually  while  the 
selected  yaw  is  set  by  the  stepper  motor.  Be  sure  to  check 
that  the  tunnel  velocity  is  set  to  the  desired  speed  prior 
to  collecting  the  data.  The  latter  is  stored  by  the  program 
in  two  files,  YPCAL1.M  and  YPCAL2.M,  which  are  later  joined 
by  MATLAB  for  use  in  curve  fitting. 

B.  DATA  COLLECTION 

The  data  collection  routine  asks  for  the  probe  position 
in  X,Y,Z  format  and  stores  the  response  as  the  first  entry 
in  a  file  named  "DATA1".  Enter  the  number  of  samples  desired 
to  be  collected  at  the  data  point,  in  response  to  the 
prompt.  The  program  then  places  the  specified  number  of 
recorded  voltages  in  the  same  data  file  as  the  point 
coordinates . 
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The  program  assumes  that  the  Probe  is  aligned  with  the 
flow  and  therefore  does  not  request  yaw  and  pitch  settings, 
if  another  orientation  of  the  probe  is  desired,  it  must  be 
recorded  at  the  time  of  testing.  Yaw  is  set  by  selecting 
Option  5,  "Yaw  Probe",  from  the  main  menu. 

C.  DATA  REDUCTION 

Data  reduction  is  accomplished  by  selecting  Option  4 
from  the  CALIBRATE  opening  menu  or  by  typing  "REDUCE"  at  the 
DOS  prompt.  The  program  will  ask  for  the  name  of  the  data 
file  to  reduce  and  for  the  name  of  file  in  which  to  store 
the  results.  Program  operation  proceeds  automatically  and 
upon  completion  of  the  reduction,  the  program  asks  if  the 
user  wishes  to  perform  another  reduction.  Data  from  thiB 
routine  is  stored  with  velocity  magnitude  in  column  1,  pitch 
angle  in  column  2  and  yaw  angle  in  column  3. 

D.  GENERAL  NOTES  ON  THE  SOFTWARE 

Raw  voltage  data  was  collected  from  the  probe  using  a 
compiled  basic  program  driver  for  a  machine  language 
subroutine  supplied  with  the  A/D  Board.  The  source  listing 
of  the  Basic  Language  program  can  be  found  in  the  appendix. 


During  this  preliminary  concept  exploration  phase,  the  DMA 


capabilities  of  the  hardware  and  software  were  not  utilized 
due  to  problems  encountered  with  the  software  supplied  with 
the  A/D  board.  Specifically,  when  Microsoft  Quick  Basic 
dynamically  dimensions  arrays,  it  places  them  outside  of  the 
current  memory  segment.  The  DASH16  machine  language  driver 
supplied  with  the  A/D  board  expects  to  find  the  arrays  in 
the  current  segment.  When  this  condition  occurs,  the  system 
crashes.  It  should  be  noted  that  data  collection  rates  of  up 
to  60  KHz  may  be  achieved  by  utilizing  these  capabilities. 

As  previously  mentioned,  a  commercially  available 
program,  Matlab,  was  used  to  perform  the  least  squares  curve 
fitting  of  the  calibration  data.  Matlab  is  a  matrix 
manipulation  program  and  provided  very  rapid  reduction  of 
large  amounts  of  data. 


time  the  best  selection  of  yaw  and  pitch  angles  for 
calibration  can  be  determined. 


I 


VI.  RECOMMENDATIONS 

A.  EQUIPMENT 

THe  current  configuration  of  the  probe  positioner  / 
holder  causes  considerable  difficultly  in  aligning  the  probe 
with  the  center  line  of  the  tunnel.  This  could  be  rectified 
by  cutting  the  probe  holder  shaft  below  the  sliding  block  on 
the  positioner  and  coupling  it  together  with  a  sleeve  and 
set  screw  arrangement.  This  would  allow  for  rotation  of  the 
probe  independent  of  the  stepper  motor. 

The  present  calibration  system  requires  manual  settings 
of  multiple  pitch  angles  and  constant  monitoring  of  the 
tunnel  velocity.  The  micromanometer  used  to  determine  tunnel 

\ 

( 

velocity  also  requires  constant  attention  to  insure  drift 

t 

from  the  zero  reference  point  has  not  occurred.  This 

f 

! 

combination  of  shortcomings  increases  the  required  j 

calibration  time.  Considering  the  frequency  with  which 

i 

J 

recalibration  is  required  to  maintain  accuracy,  replacement 
of  the  present  system  with  a  DISA  Calibrator  Unit  would 

i 

i 

reduce  the  tedium  of  the  process  as  well  as  improve  j 

i 

i 

calibration  results.  I 
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B .  PROBE 


Probe  geometry  has  been  assumed  to  be  perfectly 
orthogonal  during  the  development  of  the  algorithms. 
Improved  accuracies  may  well  be  obtained  by  direct 
measurement  of  the  angles  of  the  wires.  It  has  recently 
come  to  the  attention  of  the  author  that  an  optical 
comparator  is  available  at  NPS  to  perform  these 
measurements . 

Reference  6  provides  an  indepth  investigation  of  the 
effects  of  roll  angles  on  the  accuracy  of  data  obtained  from 
the  particular  probe  being  used.  This  information  should  be 
utilized  when  actual  measurements  are  performed. 


APPENDIX 
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’*  * 

’ *  Probe  Calibration  Routine  * 

’*  Gary  J.  Selman  Ver.  1.0  * 

’*  Naval  Postgraduate  School  14  November  1987  * 

> ********************************************************************** 

DIM  DI0X(4) 

COMDN  SHARED  DIOX< )  ,LX,UXfTEMP  ’Make  available  to  all  routines 

*********************************************************************** 
’*  * 

’*  START  OF  FUNCTION  DEFINITIONS  * 

’*  * 

>********************************************************************** 


DEF  FNVelocity(Delta.P) 

STATIC  delp 

delp=- . 00202  +  . 93*Delta.P  ’Converts  delta  p  from  static 

’port  to  equivelent  pitot  static 
’delta  p 

FNVelocity  =(2.9  *  SGR<  (T01P+459.O)*delp)  )*.3048 

END  DEF 


DEF  FNDegRad ( DEG ) =DEG  *  (Pi/180) 

DEF  FNVelocityParam(El,E2fE3)=SqR<Er2+E2~2+E3“2) 
DEF  FNPitchParam(El  ,E2,E3)=E1/(SQR(E2~2+E3'S2) ) 
DEF  FNYawParam( E2 , E3 ) =E2/E3 


’Degrees  To  Radians 
’Velocity  Parameter 
’Pitch  Parameter 
’Yaw  Parameter 


********************************************************************** 

*  SUBROUTINE  SAMPLE 

*  Samples  each  wire  specified  LX  and  UX.  Note  that  these  two  values 

*  have  been  declared  common  and  are  available  to  all  routines 

*  ARRAY  :  Three  dininsion  Array  that  returns  the  mean 

*  Voltage  of  each  wire  sampled 

*  NoSamples  :  Specifies  how  many  samples  will  be  used  to  calculte 

*  the  mean  Voltage.  If  =  1,  returns  only  an  instantanious 

*  voltage 

*  FirstDim  :  Specifies  first  diminsion  of  ARRAY. 

*  DIOX(O)  :  Returned  by  DASH 16  =  Voltage  lead  by  A/D 

*  DI0X(1)  :  Returned  by  DASH16  =  A/D  Channel  read 

*  LX  &  UX  :  Lower  and  upper  channel  scan  limits 


SUB  Sample ( NoSamplesX , FirstDim , ARRAY ( 3 ) )  STATIC  ’Note  that  LX  &  UX 

’sue  shared 


DIM  SunVolt#(3) 
ERASE  SunVolt# 


STATIC  MeanVoltage# , I , J 
MeanVoltage#=0 . 0 
MDX  =  3 

FOR  1=1  TO  NoSamplesX 

FOR  J  =  LX  TO  UX 
CALL  DASH16  (MDX,  DIOX(O),  FLAGX) 

IF  FLAGXOO  THEN  PRINT  "Error  #";FLAGX;"  in  mode" 
SumVolt# ( DIOX( 1 ) ) = ( ( 5 . *DIOX( 0 ) ) /2047 ) +SumVolt# ( DIOX( 1 ) ) 
MeanVoltage#=SunVolt#(DICX( 1 ) )/I 
ARRAY(FirstDim,l,DIOX(l))  =  MeanVoltage# 

NEXT  J 
NEXT  I 

END  SUB 


’XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 


SUB  INITIALIZE. AD  STATIC 

OPEN  "DASH16 . ADR"  FOR  INPUT  AS  #1 
INPUT  #1,  DIOX(O) 

CLOSE  #1 
DIOX(l)  =  2 
DIOX( 2)  =  3 
DASH16  =  0 
FLAGX  =  0 
MDX  =  0 


'initial  with  mode  0 
'get  base  I/O  address 


’interrupt  level 
'D.M.A.  level 

'call  offset  -  always  zero 
'error  variable 
'mode  0  -  initialize 


CALL  DASH16  (MDX,  DIOX(O),  FLAGX) 

IF  FLAGXOO  THEN  PRINT  "INSTALLATION  ERROR": STOP  'Halt  on  error 
END  SUB 


>XXXXXXXXXXXXXXXXXXXXX*XXXXX*XXXXXXX*XXXX*XXX*XX***********X*********** 


SUB  GET. SCAN. LIMITS ( LX, UX)  STATIC  ’  set  scan  limits 

'Scan  limits  will  default  to  0-7  (8  channel)  or  0-15  (16  channel)  if 


you 

'skip  this  step 

INPUT  "Enter  lower  channel  scan  limit  number  :  ";LX 
INPUT  "Enter  upper  channel  scan  limit  number  :  "  ;UX 
END  SUB 

’XXXXXXXXXXX**XXX*XXXX*»**XXXXXXXXXXXXXXX**XXXXX**«********XXXXXXXXXXXX 


SUB  SET. SCAN. LDHTS(LX.UX)  STATIC 
DIOX(0)=LX 
DIOX( 1 ) =UX 
MDX  =  1 

CALL  DASH 16  (MDX,  DICK(O),  FLAGX) 
IF  FLAGXOO  THEN  PRINT  "Error  #"; 


'set  lower  limit 
'set  upper  limit 
'mode  1  -  set  scan  limits 

AGX;"  in  setting  scan  limits"  :ST0P 


END  SUB 

'X********************************************************************* 

SUB  MATLAB.OUT(ARRAY(3) )  STATIC 
STATIC  I, J,K 
DIM  Wire# (3) 

Wire$(0)="NoWires=" 

Wire$( 1 )="wirel=[" 

Wire! ( 2 ) ="wire2= [ " 
wire! ( 3 ) ="wire3=  t " 
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PRINT  "  The  program  will  generate  a  file  called  VELCAL.M” 
PRINT  "  This  file  will  contain  the  data  collected  for  each” 
PRINT  "  wire  in  a  format  suitable  for  execution  by  MATLAB." 
PRINT  "  MATLAB  will  perform  a  least  squares  curve  fit  on  th 
PRINT  "  data  by  utilizing  a  program  called  DataFit.M" 

NAME$="C: \FROBCAL\VELCAL.M" 

OPEN  NAME!  FOR  OUTPUT  AS  #2 
PRINT#2 , Wire! ( 0 ), UBOUND ( ARRAY , 3 ) 

FOR  1=1  TO  UBOUND (ARRAY, 3) 

IF  I  <>1  THEN  PRINT# 2  ,"]" 

PRINT#2  ,Wire$(I), 

FOR  J=1  TO  UBOUND ( ARRAY, 1) 

FOR  K  =  1  TO  UBOUND! ARRAY, 2) 

PRINT#2, ARRAY! J.K.I), 

NEXT  K 
PRINT#2, 

NEXT  J 
NEXT  I 
PRINT#2 , " ] " 

CL0SE#2 

END  SUB 

’****3M*M***M************«***t**S«t**«»tSt***t************«**tt»tt*«* 

SUB  Velocity. Sensitivity! No. Of .Wires)  STATIC 

INPUT  "  How  many  calibration  points  will  you  enter” ; Points 
REDIM  CAL. ARRAY  (Points, 2, No. Of .Wires) 

FOR  I  =  1  TO  No.Of. Wires 
LX  =  I  :  U%  =  I 
CALL  SET.  SCAN.  LIMITS!  LX, UX) 

CLS 

PRINT  USING  "  Set  Wire  #  perpendicular  to  flow"; I 
IF  I  =  1  THEN 

PRINT  "CENTER  THE  3  DOTS  ON  TCP  OF  THE  PROBE" 

PRINT  "PRESS  ANY  KEY  WHEN  READY: 

PRINT  "«*»**#*«***#****#**  WARNING  MS*******#**********" 

PRINT  "PROBE  WILL  ROTATE  WHEN  YOU  DO" 

HEREA:  AS  =  INKEYS  :  IF  AS  =  ""  THEN  GOTO  HEREA 

Yaw  =  -45 

CALL  YawProbe ( 1 , Yaw ) 

END  IF 

IF  I  =  2  THEN 

PRINT  "CENTER  THE  SINGLE  DOT  ON  TOP  OF  THE  PROBE" 

PRINT  "PRESS  ANY  KEY  WHEN  READY: 

PRINT  "««*«*****#*»**»****  WARNING  *******«**»******»«**" 

PRINT  "PROBE  WILL  ROTATE  WHEN  YOU  DO" 

HEREB:  AS  =  INKEYS  :  IF  AS  =  ""  THEN  GOTO  HEREB 

Yaw  =  -45 

CALL  YawProbe! 1 , Yaw) 

END  IF 

IF  I  =  3  THEN 

PRINT  "CENTER  THE  2  DOTS  ON  TOP  OF  THE  PROBE" 


PRINT  "PRESS  ANY  KEY  WHEN  READY: 

PRINT  "****  WARNING  **»*" 

PRINT  "PROBE  WILL  ROTATE  WHEN  YOU  DO" 

HEREC:  AS  =  INKEYS  :  IF  AS  =  " "  THEN  GOTO  HEREC 

Yaw  =  -45 

CALL  YawProbe ( 1 , Yaw ) 

END  IF 

CLS 

PRINT  USING  "  Wire  #  Calibration  Data  ";I 

PRINT  "  VOLTAGE  VELOCITY" 

FOR  COUNT  =  1  TO  Points 

INPUT  "Delta  P  from  Mieromanometer" ,DELP 
CALL  Sample ( 1000 , COUNT , CAL . ARRAY ( ) > 

CAL . ARRAY( COUNT , 2 , I ) =FNVeloci ty ( DELP ) 

PRINT  USING  "  ##.####  ###.#### 

CAL . ARRAY ( COUNT , 1 , I ), CAL . ARRAY ( COUNT , 2 , I ) 

NEXT  COUNT 
Yaw  =  0 

CALL  YawProbe ( 1, Yaw) 

NEXT  I 

CALL  MATLAB. OUT( CAL. ARRAY ( ) ) 

END  SUB 

SUB  YawPi tchParam ( DELP , ALPHA , BETA )  STATIC 
CORRECT: 

INPUT  "  ENTER  PROBE  AXIS  PITCH  (DEGs)  :  ".ALPHA 
INPUT  "  ENTER  PROBE  AXIS  YAW  (DEGs)  :  ".BETA 
CALL  YawProbe (1, BETA) 

INPUT  "  TUNNEL  DELTA  P  :  " ,DELP 
INPUT  "  ANY  CHANGES  (Y/N)  :”;AS 
IF  AS  =  "Y"  OR  AS  =  "y"  THEN  GOTO  CORRECT 

END  SUB 

SUB  YawProbe (Mode*, Yaw)  STATIC 
DIM  Pointer(7) 

STATIC  N, Position 
IF  ModeX  =  0  THEN 
CLS 
N  =  3 

DATA  9,5,6,10,9,5,6,10 
FOR  I  =  0  TO  7 

READ  Pointer! I) 

NEXT  I 

PRINT  "  Insure  probe  is  aligned  with  tunnel  center  line' 
print  "  " 

PRINT  "  DO  NO  CONNECT  POWER  TO  STEPPER  MOTOR 
PRINT  "  " 

PRINT  "  Press  any  key  to  continue" 

HERE1 :  AS  = INKEYS:  IE  AS  =  ""  GOTO  HERE1 

MDX  =  13 

DIOX(O)  =  Pointer! N) 
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PRINT  USING  "DIOX ( 0 )  =  ##";  DIOX(O) 

CALL  DASH 1 6 ( MDX , DIOX (  0  )  , FLAGX ) 

IF  FLAGX  <>  0  THEN  PRINT  "Error  " : STOP 
Position  =  0 
’  CLS 

PRINT  "A  valid  code  now  exists  on  the" 
print  "input.  Power  may  now  be  connected" 
print 

print  "PRESS  ANY  KEY  TO  CONTINUE" 

HERE2:  A$  =  INKEYS  :  IF  AS  =  ""  GOTO  HERE2 


ELSE 

MDX  =  13 

Direction  =  0 

WHILE  (Yaw  <>  Position) 

PRINT  USING  "Yaw  =  ##.##  Position  =  ###.##" ;Yawf Position 

IF  Yaw  >  Position  THEN  Direction  =  -1  ELSE  Direction  =  1 
N  =  N  +  Direction 

IF  N  <  0  AND  Direction  =  -1  THEN  N  =  7 
IF  N  >  7  AND  Direction  =  1  THEN  N  =  0 
DIOX(O)  =  Pointer(N) 

CALL  DASH16(MDX,DIOX(0) , FLAGX) 

FOR  Delay  =  1  TO  5000 
NEXT  Delay 

Position  ^Position  -  Direction*15 
WEND 

PRINT  USING  "Yaw  =  ##.##  Position  s  ###.##"; Yaw, Position 
END  IF 

END  SUB 


' ********************************************************************** 
’*  * 

’  *  START  OF  MAIN  PROGRAM  * 

’  *  * 

’ ***************************************** ***************************** 


CALL  INITIALIZE. AD  ’  Initialize  A/D  Converter 

CALL  YawProbe(O.l) 

SCREEN  8,0,1,1:WIDTH  80:CLS:KEY  OFF 
Pi  =  3.14159263 

MAIN: 

CLS 

PRINT  "Do  you  want  to  1.)  King’s  Law  Calibrate  " 

PRINT  "  2 . )  Yaw  /  Pitch  calibrate  " 

PRINT  "  3.)  Collect  Data  " 

PRINT  "  4.)  Reduce  Data" 

PRINT  "  5 . )  Yaw  Probe 

PRINT  "  6.)  Quit  " 

INPUT  "Enter  Your  Choice"; A 

ON  A  GOTO  KING, YAW. PITCH, COLLECT, REDUCE, SetYaw, QUIT 
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INPUT  "Enter  Ambient  Temperature" ; TEMP 
INPUT  "How  Many  Wires  No. Wires 
CALL  Velocity . Sensitivity! No .Wires ) 

GOTO  MAIN 
YAW. PITCH: 

SCREEN  8,0, 2,1 

PRINT  "  PRESS  ANY  KEY  TO  CONTINUE" 

SCREEN  8,0, 1,1 
CLS 

DIM  E( 1 , 1 ,3) 

OPEN  "C:\PR0BCAL\YFCAL1.M"  FOR  OUTPUT  AS  #3 
OPEN  "C:\PROBCAL\YFCAL2.M"  FOR  OUTPUT  AS  #4 
PRINT#3, "ypl=[" , 

PRINT#4,"yp2=[", 

PRINT  "This  routine  assures  that  3-D  probe  calibration  is  in 

progress" 

INPUT  "  AMBIENT  TEMPERATURE  :  "  ,TEMP 
CALL  GET. SCAN. LIMITS! LX ,UX) 

CALL  SET. SCAN. LIMITS! LX, UX) 

AGAIN: 

CALL  YawPi tchParam ( DELP , ALPHA , BETA ) 
counts 1 . 

CALL  Sample! 1000, count, E! )) 

E1=E( 1,1,1) 

E2=E( 1,1,2) 

E3=E( 1, 1,3) 

CalVeloci ty=FNVeloci ty ( DELP ) 

X=FNVeloci tyParam! El , E2 , E3 ) 

Y=FNPitchParam(El ,E2,E3) 

Z=FNYawParam(E2 , E3 ) 

CalVelocitysFNVelocity ( DELP ) 

SCREEN  8,0, 2, 2 

PRINT  "  Velocity  Yaw  Pitch  Sensor  Voltage" 

PRINT  USING  "  ##.####  #*.*##  ##.### 

## . ####" ;CalVelocity, 

BETA, ALPHA, El 

PRINT  USING  "  ## . ####" ; E2 

PRINT  USING  "  ##.####" ;E3 

HERE:  AS=INKEY$ :  IF  A$='"’  GOTO  HERE 

SCREEN  8, 0,1,1 

ALPHA= FNDegRad! ALPHA ) 

BETA=FNDegRad  ( BETA ) 

PRINT#  3 ,  CalVelocity ,  ALPHA ,  BETA ,  X ,  Y 
PRINT#4,Z,E1,E2,E3 

INPUT  "DO  YOU  WANT  TO  ENTER  CALIBRATION  ANOTHER  POINT"  ;A$ 

IF  A$  =  "Y"  OR  A)  :  "y"  GOTO  AGAIN 
FRINT#3 , " ] " 

PRINT#4 , " ] " 

CLOSE  #3 
GOTO  MAIN 


COLLECT : 


CLS 

CALL  GET. SCAN. LIMITS ( LX, UX) 

CALL  SET. SCAN. LIMITS ( LX, UX) 

DataPoint  =  1 
CollectAgain: 

INPUT  ”  Enter  X,Y,Z  coordinates  of  sample" ,X,Y,Z 
INPUT  "  Enter  Number  Of  Samples  :  ".COUNT 
FILE!  =  "C:\PROBCAL\DATA"+STR$( DataPoint) 

OPEN  FILES  FOR  OUTPUT  AS  #2 
REDIM  SDATA( COUNT ,1,3) 

REDIM  S1DATA(C0UNT,3) 

S1DATA( 0 , 1 ) =X:S1DATA( 0,2) =Y: S1DATA( 0 , 3 ) =Z 
PRINT  #2  ,X,Y,Z 
FOR  I  =  1  TO  COUNT 

CALL  Sample (1,1, SDATA( ) ) 

NEXT  I 

FOR  1=1  TO  COUNT 
FOR  J  =  1  TO  3 

S1DATA( I , J ) =SDATA( I , 1 , J ) 
PRINT#2,S1DATA(I,J), 

NEXT  J 
PRINT#2 , 

NEXT  I 

CLOSE  #2 

DataPoint  =  DataPoint  +  1 

PRINT  "Do  you  want  to  :  1)  Reduce  this  data  now" 

PRINT  "  2)  Sample  another  position 

PRINT  "  3)  Quit" 

INPUT  "Enter  your  choice  :  ";A 

ON  A  GOTO  REDUCE, CollectAgain, MAIN 

REDUCE: 

RUN  "REDUCE" 

SetYaw: 

INPUT  "Enter  desired  yaw  angle  :  ",Yaw 
CALL  YawProbe ( 1 , Yaw ) 

GOTO  MAIN 

QUIT: 


't*t*t*ttttttt**tt**t**t**t***tttttt**t**tttt**t*ttttt*tt*t****tttt*t** 
’*  * 
’*  3-Wire  Probe  Data  Reduction  Routine  * 
’  *  Gary  J .  Selman  Ver  .1.0  * 
’*  Naval  Postgraduate  School  14  November  1987  * 
>  tt**t*****tt**ttt*tt*t*ttt*t***ttt**t**t*t*****tt*tttt*ttttt*ttt*t tttt 


’  *  * 

’*  START  OF  FUNCTION  DEFINITIONS  * 

’*  * 

OOMION  SHARED  Pi 

DEF  FNVelocity(Delta.P) 

STATIC  delp 

delp=-. 00202  +  .93*Delta.P  ’Converts  delta  p  from  static 

’port  to  equivelent  pitot  static 
’delta  p 

FNVelocity  =(2.9  *  S®(  (TEMP+459.0)*delp)  )*.3048 
END  DEF 


’Degrees  To  Radians 
’Velocity  Parameter 
’Pitch  Parameter 
’Yaw  Parameter 

DEF  FNacos ( X ) 

IF  X  =  0  THEN 

acos  =  Pi/2 

ELSE 

acos  =  ATN(SQR(1-X~2)/X) 

END  IF 

FNacos  =  acos 
END  DEF 

DEF  FNasin(X)  =  ATN(X/SQP( l-X'2) ) 


DEF  FNDegRad ( DEG ) =DEG  *  (Pi/180) 

DEF  FNVelocityParam(El , E2 ,  E3 )  =SQR(E1  ~2+E2~2+E3~2 ) 
DEF  FNPitchParam(El ,E2,E3)=E1/(SQB(E2"2+E3~2) ) 
DEF  FNYawParam(E2,E3)=E2/E3 


SUB 


SelectThetas( THETA ( 1 ) , Select 1 ,Select2 )  STATIC 
Select 1  =  1 


Select2  =  1 


DIM  R(3) 


FOR  I  =  1  TO  3 

IF  THETAd)  >  Pi/2  THEN 

R(I)  =  ABS(Pi/2-(THETA( I)-Pi/2) )/(THETA( I )-Pi/2) 
ELSE  ’Normilized 


R( I)  =  ABS(Pi/2  -  THETA(I) )/THETA(I)  ’ratios  indicating 


END  IF 

IF  R(I)  <  1  THEN  R(I)  =  1/R(I) 
NEXT  I 

IF  R( 1 )  >  R(2)  THEN  Select 1  =  2 

IF  R ( Select 1)  >  R(3)  THEN  Select 1  =  3 

R( Select 1)  =  1000 


'distance  from  0  or 
’90  degrees 

’Select  the  two 
’angles  farthest 
’from  0  or  90 
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IF  R  (  1  )  >  R ( 2 )  THEN  Select2  = 
IF  R(Select2)  >  R(3)  THEN  Select2  =  3 
END  SUB 


SUB  ThetaLim(DC(2),THETA(l),Sl,S2,S3,L,M,Z,LIM)  STATIC 
DIM  Qi<3, 3) 

DIM  FF{ 3) 

IJK  =  DC( SI , 1 ) *DC(S2 , 1 )  +  DC(S1,2)*DC(S2,2)  +  DC(S1 , 3 ) *DC(S2 , 3 ) 

LI  =  DC(S1,1)*DC(S1,2)/M 

L2  =  (L/M~2)*(DC(Sl,2r2+DC(Sl,3r2) 

LL  =  2*(L1  -  L2) 

B  =  (-4*Z*DC{S1,3)',2) 

BIB  =  (LL~2  -  4*Z* (DC(S1 , 2 ) ~2+DC( SI , 3 ) ~2 )/M~2 ) 
dim  =  SQR(B/B1B)*SIN(THETA(S1) ) 
lim  =  FNacos(clim) 

FOR  J  =  1  TO  3 

FF(J)=((DC(S1,J)*C0S(THETA(S1)))  -  DC(S2, J)*OOS(lim) )_ 
/SIN(lim) 

NEXT  J 

Tr  =  SIN(THETA(S1 ) )/SIN(lim) 

N=(DC(S2, 1 )*FF( 1 )  +  DC(S2,2)*FF(2)  +  DC(S2,3)*FF(3) )/Tr 
W1  =  (N/M)*DC(S1,1)*DC(S1,2) 

W2  =  (N*L/M“2)*(DC(S1,2)"2  +  DC(S1,3)A2) 

W  =  2*(W1  -  W2) 

X  =  (N~2/M~2)*(DC(Sl,2r2  +  DC(Slt3)~2)  -  DC(S1,3)^2 
W3  =  SQR(abs( W~2  -  4*Z*X) ) 

Qn(Sl , 1 )  =  -W/(2*Z) 

(^i(Sl,2)  =  (N  -  <4i(Sl,l)*L)/M 

Qn(Sl ,3)  =  (-®i(Sl,l)*DC(Sl,l)  -  Qi(Sl , 2 )*DC(S1 , 2) )/DC(Sl , 3) 
if  vect  >  1.1  then  stop 

SBETA  =  (-DC(S1,2)*OOS(THETA(S1))  -  Qi(Sl ,2)*SIN(TOETA(S1 ) ) 

BETA  =  FNas in (SBETA) 

SALPHA  =  DC(Sl,l)*OOS(THETA(Sl))4^i(Sl,l)*SIN(THETA(Sl)) 

S ALPHA  =  Pi/2  -  ( SALPHA/ COS ( BETA ) 

ALPHA  =  FNasin( SALPHA) 

Slim  =  -DC ( S3, 1)*OOS< ALPHA )*COS (BETA )_ 

-  DC(S3,2)*SIN(BETA)+DC(S3,3)*SIN(ALPHA)*OOS(BETA) 
lim  =  FNacos(Slim) 

END  SUB 

*********************************************************************** 

>********************************************************************** 
’*  * 

’*  START  OF  MAIN  PROGRAM  * 

’*  * 

>********************************************************************** 
WIDTH  80:CLS:KEY  OFF 
Pi  =  3.14159263 

REDUCE: 
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DIM  DC ( 3 , 3 ) 

DIM  $i<3,3) 

DIM  A(9) :DIM  V(3) 

DIM  B(3,14)  :DIM  Qeff(3) 

DIM  K(3,3):DIM  THETA(3) 

DIM  FF(3) 

OPEN  "C:\PROBCAL\APARAM.DAT"  FOR  INPUT  AS  #1 
OPEN  "C:\PROBCAL\BPARAM.DAT"  Ft*  INPUT  AS  #2 
FOR  I  =  0  TO  9  :  INPUT# 1 , A( I ) : NEXT  I 
FOR  I  =  0  TO  14 

FOR  J  =  1  TO  3 
INPUT#2,B(J,I) 

NEXT  J 
PRINT 
NEXT  I 
CLOSE  #1 ,#2 

OPEN  "C:\PROBCAL\KPARAM.DAT"  FOR  INPUT  AS  #2 

OPEN  "C : \PROBCAL\DIROOS "  FCR  INPUT  AS  #4  ’get  directional 

FOR  I  =  1  TO  3  'cosines 

FOR  J  :  1  TO  3 

INPUT#2,K(I, J)  'get  velocity  coeff 

INPUT#4,DC(J,I) 

NEXT  J 
NEXT  I 
CLOSE  #4 ,#2 

INPUT  "Enter  data  file  name  to  reduce  :  ",FILE$ 

INPUT  "Enter  data  file  name  to  store  results  :  ",FILE1$ 

OPEN  FILE1$  FOR  OUTPUT  AS  #2 
OPEN  FILES  PC*  INPUT  AS  #1 
open  "c:\probcal\theta.dat"  for  input  as  #5 
First  =  1 


BEGIN: 


WHILE  NOT  EOF(l) 


SkipFirst: 

INPUT  #1  ,V(1),V(2),V(3) 
IF  First  =  1  THEN 
First  =  2 
GOTO  SkipFirst 


’The  first  point  is  the 
’  data  point  location  so 
’  dont  reduce  it! 


END  IF 

PC*  I  =  1  TO  3  ’Calculate 

Qeff(I)  =  K(II1)+K(I,2)«V(I)+K(I,3)*V(I)"2  ’effective 


’velocities 


NEXT  I 

X=  FNVelocityParam(V( 1 ) ,V(2) ,V(3) )  ’  Velocity  Parameter 

Y=  FNPitchParam(V(l)  ,V(2)  ,V(3) )  ’  pitch  parameters 

Z=  FNYawBaraa(V(  1 )  ,V(2) )  ’  Yaw  parameter 

Q=A ( 0 ) +A( 1 ) *X+A<  2 ) *X" 2+A ( 3 ) *Y+A( 4 ) *Y~  2+A( 5 ) *Z+A( 6 ) *Z~ 2_ 
+A(7)*X*Y+A(8)*X#Z+A(9)*Y*Z  '  Velocity  Magnitude 

FOR  I  =  1  TO  3  ’  Caculate  Thetas 

F=Qeff(I)/Q 

PARTA  =  B(I,0)  +  B(I,1)*Q  ♦  B(I,2)*Q~2  ♦  B(I,3)*F_ 

+  B(I,4)*F~2  +  B(I,5)*Y  +  B(I,6)*Y~2  +  B(I,7)*Z_ 
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+  B(I,8)*Z~2 

PARTB  =  B( 1,9 ) *Q*F  +  B(I,10)*Q*Y  +  B(I,11)*Q*Z_ 

+  B(I,12)*F*Y  +  B(I,13)*F*Z  +  B(I,14)*Y*Z 
THETA ( I )  =  PARTA  +  PARTB 

NEXT  I 

CALL  SelectTbetas(THETA( ) ,S1,S2)  ’SI  k  S2  will  be  the  index 

’of  the  two  angles  farthest 
’from  0  or  90  degs 
IF  SI  >  S2  THEN  SWAP  S1,S2 

IF  SI  =  1  AND  S2  =  2  THEN  S3  =  3 

IF  SI  =  2  AND  S2  =  3  THEN  S3  =  1 

IF  SI  =  1  AND  S2  =  3  THEN  S3  =  2 

L=(DC(S2,3) *DC(S1 , 1 )/DC(Sl ,3 ) )-DC(S2, 1 ) 

M= ( DC ( S2 , 3 ) *DC ( S 1 , 2 ) /DC ( S 1 , 3 ) ) -DC ( S2 , 2 ) 

Z  =  (DC(S1, 1 )  -  DC(S1,2)*(L/M))A2  +  (DC(S1 ,3r2)*(  1+(L~2/M~2) ) 
CALL  ThetaLira(DC( ) , THETA ( ) ,S1 ,S2,S3,L,M,Z,lim) 

IF  THETA (S3)  >  lim  THEN  SIGN  =  1  ELSE  SIGN  =  -1 
FOR  J  =  1  TO  3 

FF(J)=(DC(S1,J)*C0S(THETA(S1>) 

DC ( S2 , J ) *OOS ( THETA ( S  2 ) )) /SIN (THETA (S2) ) 

NEXT  J 

Tr  =  SIN(THETA(S1 ) ) /S IN ( THETA (S2) ) 

N= ( DC ( S2 , 1 ) *FF ( 1 )  ♦  DC(S2,2)*FF(2)  +  DC(S2,3)*FF(3) )/Tr 
W1  =  (N/M)*DC(S1,1)*DC(S1,2) 

W2  =  (N*L/M~2)*(DC(S1 ,2)“2  +  DC(S1,3T2) 

W  =  2*(W1  -  W2) 

X  =  (N"2/M^2)*(DC(S1,2)A2  +  DC(S1,3)“2)  -  DC(S1,3)~2 
Qn(Sl , 1 )  =(-W  ♦  (SIGN*SQR(ABS(W~2  -  4*Z*X) ) ) )/( 2*Z) 

Qn(Sl ,2)  =  (N  -  (*i(Sl,l)*L)/M 

Qn(Sl,3)  =  ( -Qi(Sl , 1 )*DC(S1 , 1 )  -  ^i(Sl,2)*DC(Sl,2))/DC(Sl,3) 
SBETA  =  (-DC(S1I2)*C0S(THETA(S1))  -  Oj(S1,2)*SIN(THETA(S1) ) 
BETA  s  FNasin( SBETA) 

S ALPHA  =  DC(S1 , 1 )*006(THETA(S1 ) )+^i(Sl , 1 ) *SIN(THETA(S1 ) ) 
SALPHA  =  Pi/2  -  ( SALPHA/COS ( BETA ) 

ALPHA  =  FNasin( SALPHA) 

ALPHA  =  ALPHA* (180/Pi) 

BETA  =  BETA* (180/Pi) 

PRINT  Q, ALPHA, BET  A 
FRINT#2,  Q, ALPHA, BETA 
input  "hit  a  key";A$ 

WEND 
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